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o 

The decay B s — > fx + fi is one of the milestones of the flavor program at the LHC. We 
reappraise its Standard Model prediction. First, by analyzing the theoretical rate in 
the light of its main parametric dependence, we highlight the importance of a complete 
evaluation of higher-order electroweak corrections, at present known only in the large-mt 
limit, and leaving sizable dependence on the definition of electroweak parameters. Using 
insights from a complete calculation of such corrections for K — > 'kvv decays, we find 
a scheme in which NLO electroweak corrections are likely to be negligible. Second, we 
address the issue of the correspondence between the initial and the final state detected 
by the experiments, and those used in the theoretical prediction. Particular attention is 
devoted to the effect of the soft radiation, that has not been discussed for this mode in 
the previous literature, and that can lead to 0(10%) corrections to the decay rate. The 
"non-radiative" branching ratio (that is equivalent to the branching ratio fully inclusive of 
bremsstrahlung radiation) is estimated to be (3.23 ±0.27) x 10~ 9 for the flavor eigenstate, 
with the main uncertainty resulting from the value of fs a , followed by the uncertainty due 
ON to higher order electroweak corrections. Applying the same strategy to Bj — > we 

find for its non-radiative branching ratio (1.07 ± 0.10) x 10 -10 . 

OO 
O 

1 Introduction 

The rare decay B s — > provides one of the best probes of the mechanism of quark- 

flavor mixing. Within the Standard Model (SM) this transition is mediated by a flavor- 
changing neutral current (FCNC) amplitude, is helicity suppressed, and is characterized 
by a purely leptonic final state. The first two features amount to a double suppression 
mechanism, responsible for the extremely rare nature of this decay. The third feature 
causes it to be theoretically very clean at the same time. All these considerations make 
the rare decay B s — > a formidable probe of physics beyond the SM, especially of 

models with a non-standard Higgs sector. The B s — > decay has not been observed 

yet, but the LHC experiments are rapidly approaching the sensitivity to observe it [l]-[3] 
(see also Ref. [4|), if it occurs at the SM rate. Indeed, the present 95% C.L. bounds 
read fll[5l 



Oh 



B(B S -> < 4.1 x 10~ 9 , B{B d -> fi + fJ,~) < 8.2 x 10~ 10 . (1) 
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In view of a precise measurement of B(B S — > in the near future, it is of utmost 

importance to assess its SM prediction to the best of our knowledge. Analogous comments 
apply to Bd — > /U + /i~. 

In order to obtain a precise prediction for the B s — > [i + /i~ rate within the SM it is 
necessary both to compute the corresponding electroweak amplitude with high accuracy 
and also to assess the correspondence between the initial and the final state detected by 
experiment, and those used in the theoretical prediction. More precisely, we can identify 
three main steps in the comparison between data and theory: 

• The evaluation of the "non-radiative" branching fraction (B^). This theoretical 
quantity is the branching fraction evaluated in the absence of soft-photon corrections. 
Thanks to the results of Refs. B^ is known in the SM including next-to- 

leading QCD corrections. As a result, it is anticipated that this quantity can be 
computed with an excellent precision, up to the parametric uncertainties from the 
-B s -meson decay constant (/b 3 ), the CKM factor (|V^Vj s |), the i? s -meson lifetime 
(tb s ), and the top-quark mass (Mf), in order of decreasing impact on the branching- 
fraction error. Nonetheless, as pointed out in |10| , at the leading order in electroweak 
corrections B^ suffers from a sizable dependence on the renormalization scheme 
for electroweak parameters like sin 2 9\y among others. While two-loop electroweak 



corrections are available in the large-mj limit thanks to Ref. 10 itself, we argue 
that residual uncertainties due to the large-mt approximation are not negligible with 
respect to the level of accuracy now required for the theory prediction. We reassess 
such corrections, and propose a scheme where they are likely to be negligible. 

• The treatment of the soft-photon radiation. In the full theory (a em 7^ 0) photon 
emission inevitably occurs in this process (strictly speaking, the width of the non- 
radiative mode vanishes). The simplest infrared-safe observable is 

£ phys (£ max ) = B(B S + n 7 )| E £ 7 <£ max (2) 

namely the branching fraction including an arbitrary number of undetected photons 
with total energy in the meson rest frame less or equal to E max . As we discuss in 
section [3] (see Refs. 11 13]), 

B ph y s (E max )=u J (E max )xB^ , (3) 

where w(£ max ) is a correction factor that we can compute with good accuracy for 
£ max <C itlb s , and which is independent of possible new-physics contributions affect- 
ing B^°\ In the limit where we consider bremsstrahlung radiation only, u(E max ) is 
known with good accuracy for any value of E max , and in this limit u{E raax ) — > 1 
when E max approaches its kinematical end point (£ max — > mB s /2). The theoreti- 
cal quantity B^ can thus be identified with the branching ratio fully inclusive of 
bremsstrahlung radiation. 

• The time dependence and initial-state tagging. Since the B s is not a mass eigenstate, 
also the nature of the initial state and how the measurement is performed in time need 
to be specified. The simplest observable accessible at hadron colliders is the flavor- 
averaged time-integrated distribution. As recently pointed out in Refs. 14,15 , B s — 
B s oscillation effects do not cancel out in this quantity because of the non- vanishing 
width difference between the two mass eigenstates. This leads to a correction factor 
with respect to the decay rate computed at initial time (t = 0) that, in principle, 
may be affected by new-physics contributions. 

The first two points apply also to the Bj — > decay. On the other hand, in the Bj 

case the complication related to the last point is absent, due to the smallness of AI^. 
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In sees. [2] and [3] we proceed with a detailed discussion of the first two points in the 



case of the B s — > + decay. Results are then summarized in sec. 3.2, where they are 
combined with the third point in order to obtain the SM prediction to be compared with 
the data. The case of Bj — > fi + fi~ is presented in sec. [I| The final section consists of the 
list of the main results of our paper and the outlook for the future. 

2 The non-radiative branching ratio 
2.1 Preliminaries 

The SM expression for the branching ratio of the non-radiative decay B s — > can be 

written as (see e.g. Ref. [To]) 



»(o) _ G F 
°s,SM - n 



a em (M z ) 
An sin 2 W 



i~B a fl a m Bs mlJl ^\V t lV ts \ 2 Y 2 (x tw ,x ht ;a s ) , (4) 

m B 3 



where Y is an appropriate loop function^ which consists of Z-penguin and box-diagram 
contributions, including QCD corrections as well as the leading electroweak corrections. 
In the absence of such corrections, Y reduces to the Inami-Lim function [17] 

/ \ X ( ^ - X 3X 1 \ fr\ 

Y ^ X) = 8{— x + jr^w lnx ) > (5) 

whose argument x can be identified, in the present discussion, with 

x tw = (6) 
lvl w 

Here mj(/i) is the top-quark mass renormalized (as far as QCD corrections are concerned) 
in the MS scheme at the scale \x. 

Leaving aside the dominant parametric uncertainty due to fs 3 for the time being, two 
evident uncertainties are present in B^g M , if Y is approximated by Yq: 

• The choice of the scale fi, which is usually chosen to be 0(mt), but could be as low as 
Mw or as high as 2mj, introducing sizable uncertainty in the branching ratio. This 
unphysical dependence has been basically removed through the NLO QCD corrections 
calculated in Refs. [6]-[9]. As these corrections have been discussed at length in the 
literature we will not elaborate on them unless necessary. This chapter is closed, at 
least for this decade. 

• The choice of renormalization scheme, or equivalently of the definition of electroweak 



parameters, pointed out in 10 . As this chapter is not yet closed and this uncertainty 
has not been discussed recently in the case of the B s — > fi + n~ decay, let us have a 
closer look at this dependence. 

In order to see that this dependence is sizable, let us compare two definitions of sin 2 #yK) 
respectively in the MS and in the on-shell scheme, in which sin 2 6w is very precisely 
known [l8] 

sin 2 W (M Z ) = 0.23116(3), [sin 2 6 W ]° S = 1 - M^/M 2 Z = 0.22290 . (7) 



1 Note that the presence of a cm in the normalization of eq. Q is fictitious: we can eliminate it expressing a em / sin 2 
in terms of Gf and Mw, thereby obtaining an expression that is well defined in the limit a cm — > 0. 
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We observe that the second choice, with all other parameters fixed, implies <8^°g M by 7% 
higher than the first choice. This corresponds to a shift in the rate by 0.22 x 10~ 9 , and 
is equivalent to a shift in fg s by 8 MeV, that is larger by almost a factor of two than 
the error of the most accurate determination of this weak decay constant. Evidently 
this renormalization-scheme uncertainty has to be removed in the era of precision flavor 
physics. 

The authors of [To] made the first step in this direction by providing the result for 
two- loop electroweak corrections to Y in the large-mj limit. At the level of the branching 
ratio, this reduced the uncertainty by roughly 30%, but a warning has been made that 
this estimate could be inaccurate and the inclusion of all NLO electroweak corrections will 
be necessary when the branching ratio in question will be precisely measured. That this 
estimate could indeed be inaccurate can be seen by keeping only the leading term in mt 
in Yq, namely x/8. This estimate misses the true value of Yq by almost a factor of two. 



In what follows we will briefly summarize the findings of 10 as far as scheme dependence 
due to missing NLO electroweak corrections is concerned, using the most recent set of input 
parameters. Subsequently we will provide a preliminary solution to this problem by using 
insights from the complete NLO calculation of electroweak corrections to K — > ttvv, that 
involved the loop function Xq(x). After this calculation the remaining uncertainty in 
K — > ttvu related to electroweak effects is far below 1% and one should hope that one day 
this will also be the case for B s — > 

2.2 Renormalization-scheme dependence 

As already mentioned, we are concerned here with the dependence upon the choice of the 
renormalization scheme for electroweak corrections. We consider four different renormal- 
ization schemes which can be distinguished by the manner sin 2 9w and the top-quark mass 
are renormalized. These are: 

• Two schemes for sin 2 9\y that in the formulae below will be distinguished by the 
parameter r s = 0,1: 

sin 2 MM z ): (r, = 0), [sin 2 %]° S : {r s = 1), (8) 

with their numerical values given in eq. ([7]). 

• Two schemes for the top-quark mass, distinguished by the parameter r% = 0, 1: 

m t ee m t (m t ) m ^ CD : (r t = 0) fn t ee m t (m t ) m ^ CD+EW : (r t = 1), (9) 
related via [To] 

mf = m 2 (l + &A t ( M ,x /lt )). (10) 

In the case of mt, only QCD corrections are MS-renormalized, whereas the mass is on- 
shell as far as electroweak corrections are concerned. In the case of mt, both QCD and 
electroweak corrections are MS-renormalized. We determine the QCD MS top-quark mass 
from the pole mass in Table [l] using RunDec [l9]|^] The explicit expression for A t (fi,Xht) 



can be found in Ref. 10 and has been calculated in [20] . 
We also define 

mj _ m 2 t M 2 G F m\ 



For the central value of M t in Table[l]we obtain m t (m t ) MS ' QCD = 163.2 GeV and m,(m t ) MS ' QCD+EW = 164.5 GeV. 



4 



21 



22 



G F = 1.16638 x 1(T 5 GeY~ z m Bs = 5.36677 GeV 

227(8) MeV 
1.466(31) ps 
| = 0.0405(8) 
= 5.27958 GeV 
190(8) MeV 
1.519(7) ps 
0.0087(2) 



a-£(M z ) = 127.937 
a s {M z ) = 0.1184(7) 
M w = 80.385 GeV 
M z = 91.1876 GeV 
M t = 173.2(0.9) GeV gig 
M h = 125 GeV W\ 
m„ = 105.6584 MeV 



/b s = 

r B 3 = 

m Bd = 
fs d = 
T B d = 



see text] 
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see text] 
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Table 1: Input parameters used in the determination of £>^°g M and £>^° SM 



explicit reference are taken from Ref. 18 



Quantities without an 
We do not show the errors for quantities whose 
uncertainty has a negligible impact on our branching-ratio determinations. The central value 
of fs B d corresponds to the central value of the lattice averages presented in Ref. 
the error is our estimate of the present uncertainty (see text for details). 
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while 



Concerning other parameters in eq. Q^we use the Fermi coupling Gf as determined 
from muon decay; a em (Mz) denotes the MS QED coupling renormalized at Mz; M\y,z 
are the on-shell masses of the electroweak gauge bosons. All the relevant parametric input 
is collected in Table [U 

Each of the four renormalization schemes in question is characterized by the pair (r s , rj). 
Once this pair is fixed, we know uniquely which of the parameters listed above is to be 
employed in the calculation of Y in ^ and which value of sin 2 9\y is to be used in the 
prefactor in this equation. Therefore in presenting a general formula for the function Y 
valid in all these renormalization schemes in the large mi-limit, we can trade the mass 
variables for the pair (r s ,rj). 

With this notation the loop function Y, including complete NLO QCD corrections [6j-|9] 
and two- loop electroweak corrections in the large- mj limit [To] is given in the (r s , rt) scheme 
as follows: 

Y{r s: r t ;cc s ) = Y eS {r s ,r t ) + ^4^-Y 1 {x tw ), (12) 



where 



Y cS (r s ,r t ) = Y (x {r t ))+Z t ^ (rf \x ht ) + 3 - 3r s ^J^ - r t A t (^x ht )) (13) 

8 \ sm z 6w J 

with 

x (r t ) = x t w + r t (x t w ~ x t w) (14) 
is the effective Inami-Lim function for the (r s ,rt) scheme. This expression generalizes the 

(2) 

that applied only to specific schemes. The explicit expression for r fe (xht) 
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formulae in 

can be found in Ref. 10 and has been calculated in [29] . Finally, the function Y\ , encoding 
the NLO QCD corrections, can be found in [9j. Note that Y\ is always evaluated in the 
MS-QCD scheme, that is using xtw-, whereas Yq is evaluated using x t w or xtw depending 
on the presence or not of the — Af term in Y e R . 



In the case of complete NLO electroweak corrections, the r s -dependence in eq. (13) 



would cancel, up to NNLO effects, the one of sin 9w in the prefactor in eq. rt4h. The 



corresponding rt dependence in the correction term in (12) would in turn cancel the one 
present in the leading term Yo. As evident from our formulae, where NLO electroweak 
corrections are only in the large-mt limit, this cancellation is only partial, implying left- 
over scheme uncertainties. 
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(r 8 ,r t ) 


£Sm[*10- 9 ] 


sin 


e w MS, 


m t 


OS 


(0,0) 


3.28 


sin 


e w Ms, 


m t 


MS 


(0,1) 


3.31 


. 2 

sin 


Ow OS, 


m t 


OS 


(1,0) 


3.42 


. 2 

sin 


OS, 


m t 


MS 


(1,1) 


3.45 



Table 2: Dependence of the B s g M prediction upon the choice of the renormalization scheme (r s 
for electroweak corrections as defined in the text. 



Using the central input values in Table [l] and for sin 2 9\y in eq. ^ we obtain the centra 
values for i3^°g M in the four renormalization schemes in question, that we collect in Table |2 
The central value in either of the cases has been obtained setting the QCD renormalization 
scale to fi = mt(rrit) MS,( ^ CI) . We will return to parametric uncertainties at the end of this 
Section. 

The following observations can be made on the basis of this table. 

• The main remaining uncertainty is due to the choice of the scheme for sin 2 9w- As 
already found in |10| , the inclusion of the NLO electroweak corrections in the large-mi 
limit reduced this scheme dependence from 7% to 5%, but the left-over uncertainty 
is disturbing. 

• The left-over uncertainty due to the choice of the scheme for the top-quark mass has 
been reduced to 0.9%. 

In summary the inclusion of the NLO electroweak corrections in the large-m^ limit 
reduced various scheme dependences but the left-over uncertainties are unsatisfactory. It 
is also possible that in other schemes the differences could be even larger. Finally, one 
cannot exclude the possibility that, after the inclusion of all NLO electroweak corrections, 
the removal of scheme dependence in the branching ratio would also shift significantly its 
value relatively to the two schemes considered. However, our analysis below indicates that 
for the (0, 0) scheme this appears not to be the case. 

For completeness, we mention that our results in table [2] do not include log-enhanced 
QED effects in the RG evolution of the Wilson coefficients [30]. These corrections, which 
are part of the full NLO electroweak terms, can be calculated by using the results in [3l], 



and have been included by Misiak in his estimate of -S^gM in Ref. 32 . They are found to 



affect the decay rate by about —1.4% [30j. Given the smallness of these contributions, we 
prefer not to include them in the absence of a full NLO electroweak analysis, comprising 
also the previously mentioned complete two-loop calculation of the electroweak matching 
conditions. 



2.3 Preliminary solution 

Our analysis shows that, without a complete calculation of NLO electroweak effects, only 
a very rough estimate of the scheme dependence can be made. At this stage we should 
emphasize that in all recent papers on B s — > fJ- + fJ,~ and most papers in the last decade this 
uncertainty has been omitted. This can be justified by the fact that most authors expected 
non-SM effects to modify the relevant branching ratio by a large amount, rendering any 
shift below 10% in the SM estimate irrelevant. With the recent stringent upper bound 
from LHCb, the situation changed dramatically and uncertainties of this size have to be 
taken into account. 
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Therefore, the question arises, which value for B^§ M should be quoted in the absence of 
complete NLO electroweak corrections. 

Here we would like to propose a preliminary solution to this problem. As already pointed 
out in [10] the same problem is present in K — > ttuu decays, which are theoretically even 
cleaner than B s — > /i + fi~ . These decays are governed by the Inami-Lim function Xo(xt), 
which differs from Yo(xt) only by box contributions. Therefore at large mt, where only 
the Z-penguin is relevant, the effective electroweak corrections to Inami-Lim functions are 
identical to the ones presented above. 

Now comes an important point. We are in the lucky circumstance that complete NLO 
electroweak corrections to K — > 'kvv have been calculated by Brod, Gorbahn and Stamou 
two years ago |33| . These authors considered three renormalization schemes: 

• The MS scheme for all parameters. In our terminology this is the (0,1) scheme. 

• The MS scheme for all couplings and the on-shell scheme for all masses. This is the 
(0,0) scheme. 

• The on-shell scheme for weak mixing angle and all masses and the QED coupling 
constant renormalized in the MS scheme. This is the (1,0) scheme. 

By calculating complete NLO electroweak corrections in these three schemes, they reduced 
the scheme dependence at the level of the branching ratio far below 1%, a remarkable result. 
Looking at the size of different corrections they concluded that the on-shell definition of 
masses, together with the MS definition of sin 2 9\y, our (0,0) scheme, is the best choice 
of the renormalization scheme, namely the scheme where NLO corrections are smallest 
in absolute value. Incidentally, we find that, in our B s — > case, this scheme is 

also the one that exhibits the smallest dependence, below 1%, upon the choice of the 
renormalization scale in the range [Mz,mt\. 

By inspection of their analysis for the mentioned scheme, in particular of equations (4.2) 
- (4.4) of their paper, a very simple prescription for the final result for K — > ixvv branching 
ratios (including complete NLO QCD and complete NLO electroweak corrections) emerges. 
Adapted to the B s — > /i + decay, this prescription is as follows: 

• Use eq. (4) for with 

sin 2 9 W = sin 2 6 W {M Z ) = 0.23116(3). (15) 

• Set 

Y(x t w, Xht', a s ) = Y (xtw) + Yi(x t w) = VyYo(xtw), Vy = 1.0113 , (16) 

47T 



where xtw is defined by eqs. ^ and (11). Our value of rjy agrees well with 1.012 quoted 
in |9). 

The complete electroweak corrections to B s — > will be different in the details, due 

to different box diagrams and the presence of charged leptons in the final state in place of 
neutrinos. Yet it is plausible to expect that the prescription given above could work here 
as well. 

We emphasize that the prescription described here can only be validated by a full-fledged 
NLO calculation of electroweak corrections. Indeed, in the case of K — > isvv the adherence 
of our simple prescription to the full NLO result can be checked because of the existence 



of such complete calculation, thanks to Ref. 33 . It is known from any perturbative cal 



culation, both in QCD and electroweak theory, that a particular definition of fundamental 
parameters and renormalization scale in the leading term allows to minimize NLO correc- 
tions. Here we are conjecturing, based on the argument presented at the beginning of this 
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section, that the same choice that works for K — > ttvu will plausibly work for B s — > 

as well. A complete analysis of NLO electroweak corrections for the latter decay is the 

only way to confirm our conjecture. 



2.4 Final result 

Using the prescription given above we now calculate Bf*§ M , with the plausible expectation 
(in the sense discussed above) that scheme dependence is kept to a minimum. As antic- 
ipated, in this case the dominant source of uncertainty in B^ M is of parametric nature. 
Treating all errors as Gaussian, we obtain the final result 



= (3.23 ± 0.27) x 1(T 9 , 



(17) 



which is closest to our large-mi result in the (0,0) scheme. We choose eq. (17) as our 
reference value for the non-radiative branching ratio. 

For future reference, we illustrate the impact of the various inputs via the following 
parametric expression: 



B 



(o) 

s,SM 



3.2348 x 10" 



173.2 GeV 



3.07 



fB 3 



227 MeV y V L466 P S 



VttVts 



(3.23 ± 0.15 ±0.23/ B J x 10~ 9 



4.05 x lO- 2 

(18) 



In the second line of eq. (18) we have explicitly separated the contribution to the error 



due to Jb s , which is the most relevant source of uncertainty and deserves a dedicated 
discussion. 

As pointed out in Ref. |34|, in principle one can get rid of the quadratic fs s dependence in 
B(B S — > /i + /x~) by normalizing this observable to ArriB a , thereby taking advantage of the 
relatively precise lattice results on the bag parameter of the B s —B s mixing amplitude, that 
enters the latter linearly. Moreover, this procedure removes also the dependence on the 
CKM parameters. Indeed, in 2003 this proposal reduced the uncertainty in B{B S — > fi + fi~) 
by a factor of three. However, given the recent progress in the direct determination of fs 3 



from the lattice 28,35-41 and in the determination of CKM parameters this strategy is 



no longer necessary although it gives presently a very similar result |42|. 

As far as the direct lattice determination of Jb s is concerned, an impressive progress has 



been made in the last years [28 , 35 - 41] . These results are summarized in [28] and included 
in the world average fs s = (227.6 ±5.0) MeV [43] . Using this result at face value we would 
get a total error on ^°g M of ±0.2 x 10~ 9 in eq. ( Jl7| . However, given that this average is 
largely dominated by a single determination 36 1 , and given that all the other unquenched 
estimates of fs a have errors of about ±10 GeV, we believe that a ±8 MeV error on fs s - 
that we deduce from the spread of the central values - is a more conservative estimate of 
the present uncertainty. 



3 Soft-photon corrections and the experimental branching ratio 

3.1 Soft-photon corrections 

As anticipated in the introduction, switching on electromagnetic interactions the B s — > 
jj,^ ji~ transition is unavoidably accompanied by real photon emission. On general grounds 
we can distinguish two types of radiation: bremsstrahlung and direct emission. The former 
is largely dominant for sufficiently small photon energies, can be summed to all orders in 
the soft-photon approximation, and leads to a multiplicative correction factor with respect 
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to the non-radiative rate. On the contrary, the direct-emission component vanishes in the 
limit of small photon energies and represents a background for the extraction of short- 
distance information on the B s — > H + [J,~ amplitude. A tight cut on the invariant 
mass (m 



close to mg s allows us to treat radiative corrections in the soft-photon 
approximation and to suppress the background due to the direct-emission component. 

In the soft-photon approximation (E max <C mB s /2), the correction factor defined in 
eq. ([3]), relating the photon-inclusive rate to the theoretical non-radiative rate, can be 
expressed as [11-13 



uj(E n 
um(E n 



^IB (E ma 

2-Emax 



1 + 



a,. 



7T 



(19) 
(20) 



where a em = 1/137.036 is the fine-structure constant and 



2A 



In 



/'/' 



1 + 8 



/</< 



1-/3 



ft 



Amf, 



(m 



1/2 



(21) 



The term wiB^max) takes into account the emission of an arbitrary number of real photons, 
with maximal energy in the meson rest frame less or equal to 



E„ 



rrii 



m 



2m Bs 



(22) 



together with the corresponding virtual corrections: infrared divergences of real and virtual 
contributions cancel out leading to this universal correction factor. The 0(a em /ir) term 



in eq. ( 19 ) represents the subleading model-dependent contribution due to infrared-finite 



virtual corrections and due to the residual contribution of real emission, that vanishes in 
the limit of vanishing photon energy. For E m3iX ~ 60 MeV, the universal term yields 



wib(60 MeV) » 0.89 , 



(23) 
x/tt) term 



amounting to a ~ 11% suppression of the non-radiative rate, whereas the 0(a e 
is expected to be below the 1% level, as discussed below. 

The normalization of E max in wrB(-E'max) is, in principle, arbitrary: different values lead 
to a redefinition of the 0(a em /7r) finite term in eq. (19). Following Ref. 13 , we normalize 
£ max to its kinematical limit (niB s /2) in order to minimize the residual finite corrections. 
The latter can be decomposed into the following three parts. 

I. The residual real and virtual corrections in the absence of direct couplings of the me- 
son to the photon. With the normalization adopted for uiiB(E max ), these corrections 
amount to 5a em / (Air) « 0.3%, corresponding to the electromagnetic corrections for 
the decay of a point-like meson fully inclusive of bremsstrahlung radiation [44] . 

II. The virtual structure-dependent terms (due to effective non-minimal couplings of the 
meson to the photon). These terms are model dependent; however, they must respect 
the helicity suppression of the non-radiative amplitude and do not contain large logs. 
As a result, they are expected to be of the same size as those in point I. 

III. The real contribution of the direct-emission amplitude. Since the direct-emission 
amplitude for B s — > /x + /i~7 is not helicity suppressed, it may represent a significant 
contribution if the E m3iX cut is not tight enough. However, according to Low's theo- 



rem 45 , the interference of bremsstrahlung and direct-emission amplitudes leads to 
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a correction to the rate that vanishes at least quadratically with the photon energy 
cut. From a naive dimensional analysis, the relative direct-emission contamination, 
for a given i? max cut, is 



<5de < 26 



2E„ 



a B{B s 

7T 



fjt+H 7) 



DE 



B 



(0) 
s,SM 



1/2 



(24) 



where B{B S — > /U + /^~7)de represents the genuine direct-emission branching fraction. 
According to the estimates in the literature (see Ref. |46| and references therein) 
the latter is 0(few x 10 -8 ). Then, if we assume B(B S — > ij. + h~j)de < 10~ 7 as a 
conservative estimate, we find that this relative correction is below 1% for E max < 
100 MeV. 



3.2 Connecting the experimental with the theoretical branching ratio 

In order to obtain a theoretical prediction for the decay rate accessible in experiments, the 
last point we need to take into account is the effect of the non-vanishing width difference 
Ar s , that has been measured recently rather precisely (47| . Following Ref. |14| , we assume 
that what is presently measured by the LHC experiments is the flavor- averaged time- 
integrated distribution, 



(B(B S f)) [t] = l - J\t' [T(B s (t') /) + T(B s (t') /)] , 



(25) 



where F(B s (t') — > f) denotes the decay distribution, as a function of the proper time (t'), 
of a B s flavor eigenstate at initial time (and correspondingly for B s ). Furthermore one 
defines 

r s = — = \ (rf + T L S ) , y s = Ts ~ Ts = 0.088 ± 0.014 , (26) 
tb s 2 21 s 

14 



H L 

with T s ' the total decay widths of the two mass eigenstates. As discussed in Ref. 
the time-integrated distribution is related to the flavor-averaged rate at t = by 

(B(B S -> f)) [t] = K f(t,y s )(B(B s -> /)) [t=0] = Kf(t,y s ) T{Bs ^ f) + ms ^ f) , (27) 

where (t,y s ) is a model- and channel-dependent correction factor. 

For the [i + [i~ final state (inclusive of bremsstrahlung radiation) the SM expression of 



the k* (t,y s ) factor is 15 
1 



Vst\ __*/ TH , ( Vat 



1 - e - t/TB " sinh ( — - e' t/TB " cosh 



t > T Bs 1 



1 - Vs 



(28) 

while the flavor-averaged branching ratio at t = is the quantity evaluated in the previous 
two sections. Putting all the ingredients together we then arrive at the following expression 

(B(B S -> /i + /i-(7)))f^ max] = Vs) x u(E max ) x , (29) 

for the quantity accessible in experiments. 
A few comments are in order: 
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The quantity which is more interesting for precise SM tests, and which can easily be 
affected by new-physics contributions, is B^\ The correction term u>(E max ) is insen- 
sitive to new physics, while n^^(t,y s ) can deviate from its SM expression only in the 
presence of new-physics models with new CP-violating phases and/or non-standard 
short-distance operators contributing to the B s — > /-i + fi~ amplitude [15|. Most im- 
portantly, the two correction terms w(£ max ) and y s ) need to be convoluted 
with the experimental efficiencies on E max and t, and, in principle, can even be de- 
termined experimentally up to their overall normalization (although an experimental 
determination of both these terms will become feasible only with a significant sample 
of B s —J- events). As a result, we encourage the experimental collaborations to 
directly provide a determination of B^°\ already corrected for these two terms. 

Since wib(tob s /2) = 1 and Kg^(t,y s ) » t/rss for t <C tb s , the theoretical quantity 
£>1°SM can be identified with the SM branching ratio of a flavor-tagged B s state at 
small times, fully inclusive of bremsstrahlung radiation only. We stress once more that 
the necessity to include the correction factor w(i? mEK ) does depend on the treatment 
of the electromagnetic radiation in the measurement. For instance, in the recent 
LHCb result (Tl, the signal is simulated fully inclusive of bremsstrahlung radiation 
and the correction term oj(E max ) (properly convoluted) is taken into account in the 
signal efficiency^] 

Finally, it is interesting to note that for the experimental choice of -E max applied by 
LHCb (J? max « 60 MeV) [I], and for t » tb s , the two correction terms in eq. (29) 
tend to compensate each other to a large extent. 



4 The B 



d 



decay 



The corresponding analysis of the — > ^ \i~ decay is a straightforward generalization of 
the one just presented for B s — > As far as the three items listed in the Introduction 

are concerned, the following comments suffice: 

• Our analysis of short-distance NLO QCD and NLO electroweak corrections (se c. [2| 
remains unchanged. In particular our conjecture, summarized by eqs. (15) and (16), 



applies identically in this case. What is trivially modified in the basic expression in 
eq. Q are the initial-state constants iriB s , tb s , Jb s and the CKM coupling, as now 
the index s is replaced by d. These new input parameters are given in Table [TJ 

The soft-photon corrections (sec. [3]) remains likewise unchanged, as the B s and B^ 
masses are very close to each other. 

The effect of Ar^ 



see sec. 



3.2) is negligible. 



Thus the final expression in eq. ( |29| ) is replaced by 

(B(B d ^^^( 7 )))^ Einax] =u;(E n 



x B 



(o) 

d,SM 



(30) 



Furthermore, using the input in Table [T] we find as the analogues of eqs. (17) and (18) 
the following results 



B 



(o) 

(ISM 



(1.07 ±0.10) x 10 



-10 



(31) 



-10 



Mi 



173.2 GeV 



3.07 



190MeVy V L519 P S 



VttV td 



1.07 ±0.05 ±0.09 



Sb 



x 10 



-10 



1.7 x lO- 3 

(32) 



We thank Tim Gershon and Matteo Palutan for useful discussions regarding this point. 
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The result in eq. (31) is our reference value for the non-radiative branching ratio of the 
Bd —> decay. In addition, similarly as for the B s case, eq. (32) illustrates the 

impact of the various inputs on the quoted central value and error. In the second line of 
this equation we have explicitly separated out the contribution to the error due to /s d , 
which is the most relevant source of uncertainty. 

Finally, eqs. (18) and (32) translate straightforwardly into a prediction for the ratio 
of the non-radiative branching ratios, B(B S — > fi + (j,~)/B(B(t — > Using /b s / fB d = 



1.195 from the separate constants in our table, and indicating its relative error as a r f 



one easily finds 



,(0) 
s,SM 

w~ 

d,SM 



30.35 1 ± 0.06 ± 2a 



s/d 



fs/d' 



(33) 



5 Summary 

In the present paper we have presented a comprehensive discussion of all the effects that 
are expected to have a significant impact on the SM prediction of B(B S — > By 
this we mean that we expect residual uncertainties to be negligible with respect to the 
foreseeable experimental accuracy. In particular we have discussed the effects of NLO 
electroweak corrections, and of the correspondence between the theoretical branching ratio 
and the experimental one, focusing on the effect of soft bremsstrahlung photons. Our main 
messages from this analysis are as follows: 

• The main uncertainty in the prediction of the non-radiative branching ratio <8^g M , by 
definition independent of soft-photon corrections, still originates in fs s , that enters 
£>s°SM quadratically. However, the impressive progress made by lattice QCD eval- 
uations in the last two years makes this error as low as 0(5%) at the level of the 
branching ratio. 

• At this level of accuracy it becomes essential to perform a complete calculation of 
NLO electroweak corrections to <8^°g M , that in this case are at present known only 
in the large-mt limit [10]. An explicit evaluation of these corrections is the only 
means by which the renormalization-scheme dependence due to the scheme choice 
for electroweak parameters like sm 2 9w-, can be reduced to a really negligible level. 
Using the large- limit approximation we estimate the present error due to unknown 
NLO electroweak corrections to be ±3%. The recently performed complete NLO 
analysis of these corrections in the case of K — > ttvu decays reduced the corresponding 
uncertainty down to per mil level [33] . 

• Anticipating the structure of complete NLO electroweak corrections in B s — > /i + pT 
to be similar to the known case of K — > ttvv, we have conjectured that the most 
reliable value for £>^°g M can be obtained by choosing sin 2 6w in the MS scheme, the 
top mass in the MS scheme only as far as QCD corrections are concerned, and taking 
the short-distance function to be the sum of the LO one and of QCD corrections only. 
This simple prescription leads to a prediction for the i3^°g M 

B %u = B ( B s ^ + /0 SM = (3-23 ± 0.27) x 10~ 9 . (34) 

This is our reference value for the non-radiative branching ratio. This result is lower 
by 2% to 7%, depending on the scheme considered, with respect to the estimates 
including NLO two-loop electroweak corrections in the large mj-limit, implying an 
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anticipated significant role of missing NLO electroweak corrections. Nonetheless, we 
have argued that, within our prescription, electroweak corrections are plausibly tiny. 



In connection with this prediction, formula (18) should allow to monitor how the 
central value for B(B S —> changes as a function of its main parametric depen- 

dencies. 

We have addressed the issue of the correspondence between the initial and the final 
state detected by the experiments, and those used in the theoretical prediction. In 
particular, we have focused on the effect of the soft radiation, that has not been 
discussed for the B s — > fJ, + n~ mode in the previous literature, and that can lead 
to 0(10%) corrections to the decay rate. We have argued that, if the sum of the 
energies of the undetected photons is small enough, the dominant effect is due to 



the correction factor in eq. (20), and we have discussed the expected magnitude of 



residual effects. This correction may provide a useful comparison yardstick against 
a more accurate Monte Carlo estimate, where non-uniform experimental efficiencies 
are properly taken into account. 

Including the effect of B s — B s oscillation, recently pointed out in 14 15 and also of 



0(10%), we arrive finally at the relation (29), that connects the theoretical with the 
experimental branching ratio. 

A completely analogous procedure, applied to the Bj — > fi + fi~ decay, leads to 



B 



(o) 

d,SM 



B(B d ^n+fj,-) 



SM 



(1.07 ±0.10) x 10 



-10 



(35) 



allowing, along with formula (32), to monitor how the central value for B{Bd 
/U + changes as a function of its main parametric dependencies. 

We are looking forward to a precise measurements of B(B S — > and B(B, 



), with the hope that they will disagree with S^°g M 
respectively. 



and M°o M in eqs. (|34|) and (|35|), 



J d,SM 
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